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ABSTRACT 

The  operation  of  replenishment  at  sea  is  investigated  for  three 
supply  ships  and  L  combatants  using  queueing  theory  concepts  and  a 
random  walk  model  in  three  dimensions.  The  distribution  for  total 
replenishment  time  given  the  initial  number  of  combatants  to  be 
replenished  by  each  supply  ship  and  a  specified  cyclic  order  of 
replenishment  is  expressed  in  terms  of  its  Laplace  transform  under 
the  assumption  of  independent,  exponential  service  times  for  each 
supply  vessel.  A  method  for  counting  all  possible  sequences  of 
replenishment  is  not  found,  but  some  preliminary  counting  techniques 
are  developed  which  may  be  useful  in  its  eventual  determination. 
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I.   INTRODUCTION 

A  major  aspect  of  modern  naval  policy  is  the  stationing  of 
task  forces  throughout  the  world,  and  the  supply  and  replenishment 
of  these  task  forces  is  a  critical  problem  with  which  naval  plan- 
ners are  faced.  It  is  vital  for  ships  at  sea  to  be  provided  with 
a  reliable  flow  of  supplies  for  them  to  maintain  a  continuous 
readiness  posture. 

Forces  afloat  are  presently  supplied  with  fuel,  food,  ammunition, 
and  other  bulk  goods  primarily  by  large  supply  ships  which  accompany 
task  forces  during  their  operations  and  may  periodically  proceed  to 
bases  ashore  to  replenish  themselves.  Transfer  of  goods  at  sea  is 
largely  accomplished  by  combatants  coming  alongside  these  support 
ships  and  receiving  supplies  by  use  of  various  rigs.  The  nature  of 
the  primary  alongside  method  of  transfer  and  its  consequent  grouping 
of  numbers  of  ships  over  a  period  of  time  increases  the  vulnerability 
of  all  participating  vessels.  Therefore,  time  of  replenishment  of 
the  whole  combatant  task  force  is  of  primary  concern.  The  problem 
of  scheduling  a  replenishment  operation  of  a  combatant  task  force 
by  a  group  of  supply  ships  such  as  oilers,  ammunition  ships,  and 
other  support  vessels  has  been  studied  by  McCullough  [3] ,  Gordon 
and  Copes  [2] ,  Patterson  [5] ,  and  Waggoner  [6] .  Waggoner  built  a 
mathematical  model  to  study  the  scheduling  of  a  fixed  number  of 
combatants  by  two  supply  ships  and  computed  the  distribution  of  total 
replenishment  time  of  all  combatants.  It  is  the  purpose  of  this  paper 
to  consider  the  analogous  problem  with  three  supply  ships. 


II.   FORMULATION  OF  THE  PROBLEM 

In  general  the  problem  may  be  considered  from  a  queueing  theory 
approach  with  each  supply  ship  representing  a  service  facility  and 
each  combatant  a  customer.  It  may  or  may  not  be  necessary  for  each 
customer  to  visit  each  facility,  depending  upon  the  needs  of  the 
customer  and  the  number  of  duplicate  facilities.  Waggoner  [6]  and 
Milch  and  Waggoner  [4]  employ  a  random  walk  model  for  study  of  a 
two  facility  queue  which  operates  in  parallel  and  series  simultaneously 
with  a  finite  number  of  customers.  Such  a  model  is  applicable  to  the 
replenishment  situation.  It  is  the  purpose  of  this  paper  to  extend 
the  random  walk  model  to  the  case  of  three  supply  ships  under  the 
following  assumptions: 

(i)  A  finite  number  of  combatants  are  initially  present, 

and  each  must  be  served  once  by  all  three  supply  ships. 
(ii)  The  service  times  at  each  supply  ship  are  independent, 
exponential  random  variables  with  parameters  X,  y,  and 
v  for  ships  I,  II,  and  III,  respectively,  and  each  one 
serves  only  one  combatant  at  a  time, 
(iii)  Service  is  simultaneous  and  consecutive.  When  a  combatant 
has  been  replenished  in  one  queue  he  moves  to  the  next 
queue  in  cyclic  order.  Each  facility  continues  in  operation 
until  it  has  served  all  combatants, 
(iv)  No  additional  combatants  join  the  queues  from  outside 
the  system  after  the  commencement  of  operations, 
(v)  When  a  combatant  has  been  served  by  all  three  supply 
ships  he  departs  the  system. 
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(vi)  Transit  times  between  queues  are  considered  negligible. 
Figure  1  illustrates  the  initial  positions  of  all  vessels  and 
the  order  of  flow  from  one  queue  to  another. 
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Figure  1 


III.   THE  THREE  DIMENSIONAL  MODEL 


Milch  and  Waggoner  [4]  have  represented  an  analagous  two  facility 
queueing  system  by  a  random  walk  in  a  two  dimensional  lattice.  The 
random  walk  begins  at  the  origin  at  time  zero  and  proceeds  by  horizontal 
and  vertical  steps  to  the  point  (L,L)  where  L  =  M  +  N,  and  M  and  N  are 
the  initial  number  of  customers  in  queues  I  and  II,  respectively.  A 
horizontal  step  corresponds  to  completion  of  servicing  in  one  queue, 
and  a  vertical  step  is  a  completed  service  in  the  other.  Howevera 
all  paths  to  (L,L)  are  not  possible  since  if  X(t)  and  Y(t)  represent 
the  amount  of  customers  served  in  each  respective  queue  at  time  t, 
then  X(t)  <  M  +  Y(t)  and  Y(t)  <  N  +  X(t) .  Figure  2  shows  a  path  in 
such  a  representation. 

Y  * 

(L,L) 


X(t)  =  L 


N 


Figure  2 

In  the  two  dimensional  model  the  distribution  of  total  operation 
time,  T,  is  fou.Td  in  terms  of  its  Laplace  transform. 
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If  three  supply  ships  are  present  as  in  Figure  1  a  three 
dimensional  random  walk  may  be  formulated  which  begins  at  the 
origin  at  time  zero  and  proceeds  by  unit  steps  at  a  time  in  one 
of  the  three  axial  directions  to  the  point  (L,L,L)  where  L  =  M  +  N  +  P, 
and  M,  N,  P  are  the  initial  number  of  ships  in  the  first/  second,  and 
third  queue,  respectively.  All  paths  to  (L,L,L)  do  not  represent  a 
replenishment  operation,  however,  since  if  X(t) ,  Y(t),  and  Z (t)  are 
the  number  of  ships  served  in  queues  I,  II,  and  III  respectively, 
then  X(t)  <  M  +  Z  (t) ,  Y(t)  <  N  +  X(t)  ,  and  Z  (t)  <.  P  +  Y(t)  are 
constraints  that  must  hold  during  the  whole  operation.  These  con- 
straints define  boundary  planes  for  the  three  dimensional  lattice. 

The  random  walk  paths  are  further  bounded  by  X(t)  =  L,  Y(t)  =  L, 
Z (t)  =  L.  For  the  purpose  of  clarity  in  further  development  these 
planes  and  their  intersections  will  be  catagorized  as  follows: 
Boundary        Constituents 

A  x=z+M,  z=L,  and  their  intersection 

B  y=x+N,  y=L,  and  their  intersection 

C  z=y+P,z=L,  and  their  intersection 

D  Intersection  of  A  and  B 

E  Intersection  of  A  and  C 

F  Intersection  of  B  and  C 

Note:    D  is  not  regarded  as  part  of  A  or  B.  Similar  statements  may 

be  made  about  E  and  F. 
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IV.   DISTRIBUTION  OF  TOTAL  REPLENISHMENT  TIME 

A.   PATH  PROBABILITIES 

The  first  move  of  the  random  walk  is  made  one  unit  in  an  axial 
direction  with  probabilities  A/A  +  y  +  v,  y/A  +  y  +  v,  and  v/A  +  y  +  v 
for  axes  X,  Y,  and  Z,  respectively,  since  these  are  the  probabilities 
of  stations  I,  II,  or  III,  respectively,  completing  replenishment  of 
the  first  combatant  under  the  assumption  of  independent,  exponential 
service  times.  Due  to  the  memoriless  property  of  the  exponential 
distribution  these  probabilities  are  true  at  all  steps  along  the  random 
walk,  except  at  points  on  the  boundaries  whose  probabilities  are  as 
follows : 


Boundary 

Axial  Direction 

Probability 

A 

Y 

\i/\l    +  V 

Z 

v/y  +  v 

B 

X 

A/A  +  v 

Z 

v/A  +  v 

C 

X 

A/A  +  y 

Y 

y/A  +  y 

D 

Z 

one 

E 

Y 

one 

F 

X 

one 

The  total  time  for  the  replenishment  is  the  time  needed  for  the 
random  walk  to  proceed  from  the  origin  to  the  point  (L,L,L)  which  is 
the  sum  of  the  times  needed  to  make  each  of  the  3  L  steps  required  for 
the  walk.  The  distribution  of  time  for  any  of  the  steps  not  originating 
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on  a  boundary  is  the  minimum  of  three  independent,  exponential 
distributions  and  therefore  is  itself  exponential  with  parameter 
A  +  y  +  v  .  Similarly  the  duration  of  steps  originating  on  bound- 
aries A,  B,  and  C  is  exponential  with  parameters  y  +   v  ,  A  +  v,  and 
A+  y,  respectively.  The  distribution  of  steps  originating  on  bound- 
aries D,  E,  and  F  is  exponential  with  parameters  v,  y,  and  A 
respectively. 

Let  Q(i,j  ,k,SL ,m,n)  be  a  route  that  has  exactly  i,  j,  k,  £,  m,  n 
positions  on  boundaries  A,  B,  C,  D,  E,  F,  respectively.  Then 
Q(i,j  ,k,£,m,n)  has  3L-  (i  +  j+k+Jl+m  +  n)  positions  along  the 
path  that  are  not  on  a  boundary.  The  total  replenishment  time,  T, 
is  composed  of  the  following  random  variables,  given  that  Q(i,j ,k,£,m,n) 
occurs: 

T  =  A.  +B.+C  +D+E  +F  +G  (4-1) 

1    j         K    I        m    n    g 

where  A.  is  the  sum  of  i  steps  each  of  which  originates  on  boundary  A, 

and  B.,  C  ,  D„ ,  E  ,  and  F  are  similarly  defined.  G  is  the  sum  of 
j       k  I       m      n  J  g 

g  steps  none  of  which  originates  on  any  boundary  where  g  =  3  L  -  i  - 
j-k-£-m-n. 

Path  Q(i,j,k,£,m,n)  contains  L  steps  in  each  of  the  three  axial 
directions.  If  the  L  steps  in  the  x  direction  are  examined,  exactly 
n  of  them  will  begin  on  boundary  F  each  with  probability  one.  Some 
number ,  j , ,  will  begin  on  boundary  B  each  with  probability  A/A  +  v 
and  some  number ,  k, ,  begin  on  boundary  C  each  with  probability  A/A  +  y . 
The  remaining  L  -  n  -  j,  -  k,  steps  will  begin  from  a  non-boundary 
position  each  with  probability  A/A  +  y  +  v .  Similar  statements  may  be 
made  concerning  the  steps  in  the  other  two  axial  directions,  and  the 
results  are  tabulated  below. 
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Axis 
Direction 

Number 
of  Steps 

Step 
Probability 

X 

n 
kl 

one 

A/A  +  v 
A/A  +  y 

L  -  (n  +  jx  +  kx) 

A/A  +  y  +  v 

y 

m 
k2 

one 

v/v   +  v 
u/a  +  y 

L  -  (m  +  i-  +  k2) 

y/A  +  y  +  v 

z 

one 

v/y  +  V 
v/A  +  v 

l  -  (a  +  ±2  + 

J2> 

v/A  +  |j  +  v 

where  i  = 

h  +  i2'  j  =  ji  + 

V 

k  = 

kl 

+  k2.  Then 

P(O)    =r(: 

L.i  .k.£  .m.n)  (— ■ — — 

4  L" 

(rrt-j 

1+k 

i}  r  X  ijl 

A+y     A+y+v  y+v     A+y 


(   v_L-(.+i2+J2)   v   i2    v   j2 

A+y+v  y+v      A+v 

P(Q)  =  r(i,j,k,*,n.,n)  (i±H±V  (*±£V 

, A+y+v. k  ,   v   vL-ft  ,   y   .  L-m  ,   A   .  L-n 

(    A+y  '   1A+m+vj     lA+y+v''     ^A+y+v' 

where  Q  =  Q(i,j,k,£,m,n)  and  r  (i,j  ,k, it, va,n)  is  the  number  of  paths 
from  the  origin  to  the  point  (L,L,L)  which  have  i,j,k,£,m,n  points 
in  common  with  boundaries  A,  B,  C,  D,  E,  F,  respectively. 
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B.   TRANSFORM  OF  THE  DISTRIBUTION 

Fran  equation  (4-1)  the  following  Laplace  transform  can  be 
tabulated : 


Random  Variable 

A. 

1 

B. 
3 

D* 

E 
m 

F 
n 


Laplace  Transform 
(y+v/y+v+s)1 
(X+v/X+v+s) 3 
(X+y/X+y+s)k 
(v/v+s)£ 
(y/y+s) 


m 


n 


(X/x+s)' 
(X+y+v/X+y+v+s)g 


The  Laplace  transform  of  the  total  replenishment  time,  T,  is: 


f*T(s) 


-ST. 


E(e  Si)  =  I         E(e  S1|q)  p(Q) 
i,j,k 
£,m,n 


_      ,.  .  ,       >  ,  Xyv  NL  ,X+y+v   N  3L 
f_(s)  =  E    r(i,],k,l,m,n)  (,  ,M,  )  (-——r^r) 


* 

:T 


i,j,k, 
£,m,n 


'X+y+v    X+y+v+S 


,X+y+v+Svi  ,X+y+y+Sv  j  ,X+yHKH-s_.k 
Vy+v+s     lX+v+s   '   (A+y+s   ' 


(4-3) 


,X+y+v+Sv  &     ,A+y+v+s>,m    /X+y+v+s.n 


v+S 


y+s 


A+S 
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v-     COUNTING  TECHNIQUES 

A.  DISCUSSION 

An  expression  for  r (i,j ,k,£ ,m,n)  has  not  been  developed  in  this 
paper.  However,  this  chapter  demonstrates  several  counting  techniques 
for  the  two  supply  ship,  two  dimensional  random  walk  which  may  be  use- 
ful in  the  eventual  determination  of  r(i,j,k,Jl,m,n).  Additionally, 
for  the  three  dimensional  model  an  expression  is  found  for  the  number 
of  paths  proceeding  from  the  origin  to  (u,v,w)  which  touch  no  bound- 
aries. In  what  follows,  the  binomial  coefficient,  (  )  =  n!/rl (n-r) ! 

B.  TWD  DIMENSIONAL  COUNTING 

1.  The  Number  of  Paths  Not  Touching  Boundaries 

First  a  method  will  be  formulated  for  counting  the  number  of 
paths  proceeding  from  (0,0)  to  (u,v)  which  do  not  touch  or  cross  either 
boundary  y  =  x  +  a,  (a  >  0) ,  or  boundary  y  =  x  +  b,  (b  <  0) .  This 
number  is  equal  to  the  total  (   )  paths  going  to  (u,v)  less  the  number 
of  paths  which  touch  or  cross  either  boundary.  By  the  Reflection 
Principle  (see  page  70  of  [1] ) ,  the  number  touching  or  crossing  y=x  +  a 
is  (  _  )  and  the  number  touching  or  crossing  y  =  x  +  b  is  (v_k)  •  Since 
in  both  of  these  quantities  paths  touching  or  crossing  both  boundaries 
are  included,  the  number  of  such  paths  must  be  added  back.  This  quantity 
may  be  computed  by  dividing  it  into  two  mutually  exclusive  categories  — 
those  paths  which  touch  or  cross  boundary  y  =  x  +  a  first,  and  those 
paths  which  touch  or  cross  boundary  y  =  x  +  b  first.  Such  paths  will 
be  counted  by  a  double  application  of  the  Reflection  Principle.  In 
Figure  3  consider  path  1  which  falls  into  the  first  of  the  two  categories. 
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The  first  reflection  is  with  respect  to  line  y  =  x  +  a.  Both  the 

boundary  y  =  x  +  b  and  that  part  of  path  1  beyond  its  first  carman 

point  with  y  =  x  +  a  are  reflected.  The  reflected  path,  path  1' , 

ends  in  point  (v-a,  u+a) .  The  common  points  of  path  1  and  y  =  x  +  b 

are  reflected  into  common  points  between  path  l1  and  the  line 

y  =  x  +  2a-b  (this  being  the  reflection  of  line  y  =  x  +  b) .  Therefore, 

the  number  of  paths  in  the  first  category  is  equal  to  the  number  of 

paths  ending  in  (v-a,  u+a)  and  touching  or  crossing  the  line  y  =  x  +  2a-b. 

The  latter  number  is 

.v-a+u+a   v    ,u+v  .  ik  -\\ 

Wa-(2a-b);    Va+b;  l ,U 

by  a  second  application  of  the  Reflection  Principle.  Similarly  the 
number  of  paths  going  to  (u,v)  after  touching  y  =  x  +  b  first  is 
(u  v,) .  Such  a  path  is  path  2  in  Figure  3 .  Therefore  the  number 
of  paths  to  (U/V)  touching  neither  y  =  x  +  a  nor  y  =  x  +  b  is 

O  -  O  -  G> +  C-b> +  Cw  m 

By  extending  the  boundary  reflection  employed  in  the  preceding 
paragraph  it  is  possible  to  compute  the  number  of  paths  proceeding  to 
(u,v)  which  touch  or  cross  y  =  x  +  a,  a  >  0,  first,  y  =  x  +  b,  b  <  0, 
second,  and  y  =  x  +  c,  c  >  a,  third,  before  going  to  (u,v) .  In  Figure 
4  since  y  =  x  +  a  is  the  first  boundary  to  be  struck,  reflect  y  =  x  +  b 
and  y  =  x  +  c  with  respect  to  y  =  x  +  a  yielding  y  =  x  +  2a-b  and 
y  =  x  +  2a-c,  respectively.  Since  y  =  x  +  2a-b  is  the  next  boundary 
to  be  struck  by  the  reflected  path,  reflect  y  =  x  +  2a-c  with  respect 
to  y  =  x  +  2a-b,  yielding  y  =  x+2a  -  2b+c.  The  twice  reflected  path 
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2a-b 


Figure  3 
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(v-a+b-c,  u+a-b+c) 


V   /N 


2a-2b+c 


2a-b 


2a-c 


Figure  4 
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will  then  strike  y  -  A+2a  -  2b+c  and  then  proceed  to  (u-a+b,  v+a-b) . 
The  number  of  such  paths  is,  by  the  Reflection  Principle, 

,u-a+b+v+a-b       ,u+v    .  . 

V+a-b- (2a-2b+c)  '  V-a+b-c;  [b  3) 

This  technique  may  be  extended  to  multiple  reflections  of  any  number 
of  similar  boundaries. 

2.  The  Number  of  Paths  Touching  a  Boundary  d  Times 

Consider  now  a  technique  for  counting  the  number  of  paths  to 
(u,v)  which  touch  boundary  y=x+a,  (a  >  0) ,  d  times,  but  do  not 
cross  it,  and  do  not  touch  or  cross  boundary  y=x+b,  (b<0).  This 
number  is  equal  to  the  number  of  paths  proceeding  to  (u,v)  which  touch 
y  =  x  +  a  d  times  without  crossing  it  less  the  number  of  paths  touching 
y  =  x  +  a  d  times  without  crossing  it  that  do  touch  or  cross  y  =  x  +  b. 
In  [4]  Milch  and  Waggoner  developed  a  new  technique  called  the  Telescope 
Principle  for  counting  the  number  of  paths  proceeding  to  (u,v)  which 
touch  y  +  x  +  a  d  times  and  do  not  cross  it.  This  number  is  equal  to 

(u+v-d)  _  ,u+v-d}  ( 

W-r    {   u+a  ;  l3  q) 

To  compute  the  number  of  paths  to  (u,v)  which  touch  y  =  x  +  a 

d  times,  do  not  cross  it,  and  touch  or  cross  y  =  x  +  b,  it  is  necessary 

to  examine  several  disjoint  subsets  of  that  set  of  paths. 

If  the  number  of  paths  in  the  set  is  N  then 

N  =  N,  +  N2  +  N3  where 

N,  =  the  number  of  paths  among  these  paths  that  touch 

or  cross  y  =  x  +  b  before  touching  y  =  x  +  a 

N_  =  the  number  of  paths  among  them  that  touch  y  =  x  +  a 

d  times  without  crossing  it  before  touching  or  crossing 
y  =  x  +  b 
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d-1 
1SL,  =  Z    (the  number  of  paths  among  them  that  touch 
3   i=l 

y  =  x  +  a  exactly  i  times  and  do  not  cross  it  before 

touching  or  crossing  y  =  x  +  b) 

To  find  N,  consider  a  path  to  (u,v)  that  crosses  y  =  x  +  b 

first,  then  touches  y  =  x  +  a  d  times  without  crossing  it.  Such 

a  path  is  shown  in  Figure  5.  If  the  path  beyond  its  first  point  of 

contact  with  y  =  x  +  b  is  reflected  about  y  =  x  +  b,  the  reflected 

path  (dotted  line)  will  proceed  to  (v-b,  u+b)  and  touch  the  reflection 

of  y  =  x  +  a,  i.e. ,  the  line  y  =  x  +  2b-a,  d  times  without  ever 

crossing  it.  The  Telescope  Principle  is  then  applied  to  the  reflected 

path.  Therefore  the  number  of  such  paths  is  equal  to 

[number  of  paths  to  (v-b,  u+b-d)  hitting  neither  y  =  x+a 

nor  y  =  x+2b-a-d] 

-  [number  of  paths  to  (v-b,  u+b-d)  hitting  neither  y  =  x+a 

nor  y  =  x+2b  -  a-d+1] 

This  number  is  obtained  from  (5-2)  yielding  the  formula  for  N, : 

_  r  ,u+v-d  v    ,u+v-d,  i    r -u+v-d   v 
1  =  L  Wa-b-r    W-b;  J    L  V2a-b-l; 

,u+v-d  n    f,u+v-d  ,  _  .u+v-d    u  (c.  r, 

Wa-bn    uu-2a+3b'    lu-2a+3b+l;  J  l°  D) 

Next  the  number  of  paths,  N  ,  proceeding  to  (u,v)  which  touch 
y  =  x  +  a  d  times  but  do  not  cross  it  then  touch  or  cross  y  =  x  +  b 
is  counted.  Such  a  path  is  shown  in  Figure  6.  If  the  Telescope 
Principle  is  applied  to  this  path  about  boundary  y  =  x  +  a,  its 
equivalent  path  (dashed  line)  can  be  seen  to  meet  the  following 
criteria : 
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(i)  it  must  touch  but  not  cross  line 

y  =  x  +  a  +  d-1  first  and  then  must 
touch  or  cross  line  y  =  x  +  b  +  d. 
It  may  touch  y  =  x  +  a  +  d-1  again 
after  touching  y  =  x  +  b. 
(ii)  it  must  go  to  (u-d,v) . 
The  number  of  paths  so  constrained  is 


N  =  A-B-C+D  =  I("Tt  J  "  (Urt>] 
2  v+a-b-1    v+a-b 


where 


r,u+v-d   v    .u+v-d  n  ,   . 

[(u+2a-b-l}  "  (u+2a-b)]  (5"6) 


A  =  the  number  of  paths  to  (u-d,v)  which  touch  or 

cross  y  =  x  +  a  +  d-1  first,  then  touch  or  cross 

y  =  x  +  b  +  d.  From  (5-1)  ,  A  =  (U*V~?  ,)  . 
2  v+a-b-1 

B  =  the  number  of  paths  to  (u-d,v)  which  touch  or 

cross  y  =  x  +  a  +  d  first, then  touch  y  =  x  +  b  +  d. 
From  (5-1)  ,B  =  (£££)  . 

C  =  the  number  of  paths  to  (u-d,v)  which  touch  or 

cross  y  =  x  +  a  +  d-1  first,  touch  or  cross 

y  =  x  +  b  +  d  second,  then  touch  or  cross 

y  =  x  +  a  +  d.  From  (5-3) ,  C  =  (  «  , ,  ,,,)  . 
1  v-2a+b-d+l 

D  =  the  number  of  paths  to  (u-d,v)  which  touch  or 

cross  y  =  x  +  a  +  d  first,  touch  or  cross 

y  =  x  +  b+d  second,  then  touch  or  cross 

y  =  x  +  a+d  again.  D  must  be  added  since 

it  has  been  subtracted  in  both  B  and  C. 

From  (5-3) ,  D  =  (U+^k  J  • 
v-2a+b-d 
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Those  paths  which  proceed  to  (U/V)  and  touch  y  =  x  +  a  (d-1) 
times  before  touching  y  =  x  +  b,  touch  y  =  x  +  a  a  total  of  d  times, 
and  do  not  touch  y  =  x  +  a  +  1  are  counted  next.  Figure  7  shows  a 
path  of  this  type.  Applying  the  Telescope  Principle  to  this  path 
about  boundary  y  =  x  +  a  its  equivalent  path  (dashed  line)  can  be 
seen  to  be  constrained  as  follows: 

(i)  it  must  touch  boundary  y  =  x  +  a  +  d-2  first, 
touch  y  =  x  +  b  +  d-1  second,  and  then  touch 
y  =  x  +  a  +  d-1. 
(ii)  it  must  go  to  (u-d,v)  . 
(iii)  it  may  not  touch  y  =  x  +  a  +  d-1  until  it  has 
touched  y  =  x  +  b  +  d-1. 
(iv)  it  may  not  touch  boundary  y  =  x  +  a  +  d. 
The  number  of  such  paths  is  equal  to 

a-  -  B-  -  c  +  D-  =  eriwo)  -  (uirdK  ,> 

u+2a-b-2  u+2a-b-l 

_   [(u+v-d       .    _    (u+v-d  }]  (       } 

uu+2a-b-l;        lu+2a-b;j  ^  n 


where 


A1  =  the  number  of  paths  to  (u-d,v)  which  must  touch 

y  =  x+a+d-2  first,  touch  y  =  x+b+d-1  second,  then 

touch  y  =  x+a+d-1.  From  (5-3)  ,  A'  =  (U+^u  .^n) 
J  v-2a+b-d+2 

B1  =  the  number  of  paths  to  (u-d,v)  which  must  touch 
y  =  x+a+d-1  first,  touch  y  =  x+b+d-1  second,  then 
touch  y  =  x+a+d-1.  From  (5-3) ,  B'  =  (^a+b-d+l)  ' 

C  =  the  number  of  paths  to  (u-d,v)  which  must  touch 

y  =  x+a+d-2  first,  touch  y  =  x+b+d-1  second,  then 

touch  y  =  x+a+d.  From  (5-3),  C  =  (U+^"f,  ,_,,). 
J  v-2a+b-d+l 
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D'  =  the  number  of  paths  to  (u-d,v)  which  must  touch 
y=x+a+d-l  first,  touch  y=x+b+d-l  second,  then 
touch  y=x+a+d.  D*  must  be  added  since  it  has 

been  subtracted  in  both  B'  and  C  .  From  (5-3)  , 

,   (u+v-d   . 
u         V-2a+b-d;  * 

When  the  numbers  of  those  paths  proceeding  to  (u,v)  which 
touch  y=x+a  1,2,...> (d-2)  times,  respectively,  before  touching 
y=x+b,  touch  y=3<+a  a  total  of  d  times,  and  do  not  touch  y=x+a+l 
are  computed  similarly,  each  is  found  to  be  equal  to  (5-7) . 
Therefore  the  total  number  of  paths  in  N_  is 

/j  -.\  ,u+v-d   N   o/j  t\  ,u+v-d   v  ,  , ,  ,  x  ,u+v-d  N 
N_  =  (d-1)  (  ,~  ,  J  -2(d-l)  (  ,  „  ,  n)  +  (d-1)  (  ,-  .  ) 
3         u+2a-b-2         u+2a-b-l         u+2a-b 

This  may  be  written  as 

la   i  \    r  /U+v-d       .         ,u+v-d       , , 
N3  =   (d_1)    [ (u+2a-b-2}   "    (u+2a-b-l) ] 

-   (d-1}    "i^a-b-l*   "   (u+2a-b)]  (5"8) 

Now  N  =  N     +  N    +  N  ,  and  from   (5-5) ,    (5-6) ,  and   (5-8) 

.u+v-d     v         .u+v-d.         .u+v-d     ,         .u+v-d         * 
Wa-b-lJ        W-b;        V2a+3bJ        lu-2a+3b+l; 

,u+v-d     >         .u+v-dv         .j, .v    r , u+v-d       x        ,u+v-d  X1  ,c  n\ 

+   Wb-lJ   -   (v+a-b)   "   (d+1)    CCu«rt-J,?   "    (u+2a-b)]  (5~9) 

+  (d-D  [("*ri  j  -  (utri  i)] 

u+2a-b-2  u+2a-b-l 

Combining  the  results  of  (5-4)  and  (5-9) ,  the  number  of  paths 
to  (u,v)  which  touch  boundary  y=x+a  d  times,  do  not  cross  it,  and  do 
not  touch  boundary  y=x+b  is 
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r  ,u+v-cl         ,u+v-d. ,        r  ,u+v-d     v    _    ,u+v-d^ , 
lW-l;        {  u+a  ;J        uu+a-b-r        W-bj  J 

_r,u+v-d     ,    _    ,u+v-d         v-   _   ,,u+v-d     >    _    ,u+v-dn 
Uu-2a+3b;         lu-2a+3b+rJ         L  W-b-r         WbJJ  ^~1V) 

,    r,u+v-d       >         ,u+v-d  n        ,,  -v    r,u+v-d       v         ,u+v-d       n 
+d   [(u+2a-b-l)    "    (u+2a-b)]   "    (d_1)    [ (u+2a-b-2}   "    (u+2a-b-l)] 


C.   THREE  DIMENSIONAL  COUNTING 

The  number  of  paths  from  the  origin  to  some  point  (u,V/W)  of  the 

three  supply  ship  model  which  do  not  touch  any  of  the  boundaries  of 

Chapter  III  will  be  counted  next.  This  number  is  equal  to  the  total 

number  of  paths  from  the  origin  to  (u^w)  ,  irrespective  of  boundaries, 

less  the  number  of  paths  which  touch  at  least  one  of  the  boundaries. 

The  total  number  of  paths  from  the  origin  to  (U/V,w)  irrespective  of 

any  boundaries  is  the  number  of  possible  combinations  of  u  steps  in 

the  x  direction,  v  steps  in  the  y  direction  and  w  steps  in  the  z 

direction.  This  number  is  the  multinomial  coefficient,  (     )  = 

u,v 

(u+v+w}  m   (U+V+W)  1/uIvlwl 
u,v,w 

In  order  to  count  the  number  of  paths  from  (0,0,0)  to  (u,v,w)  which 

contact  boundary  planes,  it  will  be  necessary  to  extend  the  Reflection 

Principle  to  three  dimensions .  By  the  Reflection  Principle  (        ) , 

w+a,v,u—a 

(        ) ,  and  (        )  are  the  number  of  paths  that  have  at  least 
v-a,u+a,w       u,w-a,v+a  ^ 

one  common  point  with  the  planes  x=z+a,  y=x+a,  and  z=y+ar  respectively, 
among  the  (     )  paths  connecting  the  origin  and  the  point  (u,v,w) . 

U  /  v  /W 

If  N1  is  the  number  of  paths  from  the  origin  to  (u,v,w)  which  touch 
at  least  one  boundary  plane  then 

N'  =  E1  +  F'  +  G'  -  H1  -  I1  -  J'  +  K' 
where 
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E1  =  the  number  of  paths  touching  x  =  z  +  M 

F1  =  the  number  of  paths  touching  y  =  x  +  N 

G1  =  the  number  of  paths  touching  z  =  y  +  P 

H'  =  the  number  of  paths  touching  x  =  z  +  M  and  y  =  x  +  N 

I1  =  the  number  of  paths  touching  x  =  z  +  M  and  z  =  y  +  P 

J1  =  the  number  of  paths  touching  y  =  x  +  N  and  z  =  y  +  P 

K*  =  the  number  of  paths  touching  x  =  z+M,  y  =  x  +  N,  and  z  =  y  +  P. 

By  direct  application  of  the  Reflection  Principle 
,  m     ,u+v+w    ) 

,  =  ,u+v+w    ) 

,    ,u+v+w    , 
Vu,w-P,v+P; 

H'  is  the  sum  of  two  numbers;  the  number  of  paths  which  strike  x  =  z  +  M 

first  and  the  number  striking  y  =  x  +  N  first.  I'  and  J'  are  the  sum  of 

similar  terms.  The  formula  for  J'  is  explained  below  and  H1  and  I1  are 

listed  in  tabular  form.  The  reflection  of  the  plane  y  =  x  +  N  about  the 

plane  z  =  y  +  Pisx=z-  (N  +  P) .  Any  path  proceeding  to  the  point 

(u,v,w) ,  initially  touching  the  plane  z  =  y  +  P  (one  or  more  times) ,  and 

then  proceeding  to  the  plane  y  =  x  +  N,  is  reflected  beyond  its  first 

point  of  contact  with  z  =  y  +  P  about  that  plane.  The  reflected  path 

touches  or  crosses  the  reflection  of  y  =  x  +  N,  i.e.,  the  plane 

x  =  z  -  (N  +  P) ,  and  proceeds  to  the  point  [ (u-N) ,  (v-P) ,  (w  +  N  +  P) ] . 

The  number  of  such  paths  is  (  „   ._) . 

t^  u-N, v-P 

The  reflection  of  plane  z  =  y  +  P  about  y  =  x  +  Nisx=z-  (N  +  P). 
Any  path  proceeding  to  the  point  (u,v,w)  ,  initially  touching  the  plane 
y  =  x  +  N  (one  or  more  times) ,  and  then  proceeding  to  the  plane  z  =  y  +  P 
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is  reflected  beyond  its  first  point  of  contact  with  y  =  x  +  N  about 
that  plane.  The  reflected  path  touches  or  crosses  the  reflection  of 
z  =  y  +  P /  i.e.,  the  plane  x  =  z  -  (N  +  P) ,  and  proceeds  to  the  point 
[(u  -  (N  +  P)) ,  (V  +  N) ,  (w  +  P)]  .  Therefore  J'  =  £t^_p)  +  (^JJJ^p) 
H1  and  I1  are  similarly  computed  and  are  tabulated  below. 


Initial  Plane 
Second  Plane 
Reflected  Plane 
Reflected  Point 
Number  of  Paths 


Table  for  H1 
y=x+N 
X+z+M 
z=y-  (M+N) 

[u-N,v+M+N,w-M] 

,u+v+w  v 
^u-N,w-M; 


x=z+M 
y=x+N 
z=y- (M+N) 

[u+M,v+N,w-(M+N] 

,u+v+w  v 
^u+M,v+N; 


,    ,u+v+w  .    ,u+v+w  v 
Vu-N,w-M;    Vu+M,v+N; 


Initial  Plane 
Second  Plane 
Reflected  Plane 
Reflected  Point 
Number  of  Paths 


Table  for  I1 

z=y+P 

x=z+M 

x=z+M 

z=y+P 

y=x-  (M+P) 

y=x-  (M+P) 

[u+M,v-(M+P) ,v+P] 

[u+M+P,v-P,w-M] 

,U+V+W   v 

lu+M,w+P 

,U+V+W   v 

V-P,w-NJ 

,  _  >U+V+W   v    ,u+v+w   . 
^u+M,w+PJ    V-P7w-M; 

K  may  be  computed  in  similar  fashion  by  permiting  all  three  boundary 
planes  and  applying  the  Reflection  Principle  twice  after  each  permutation. 
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